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METHOD FOR RADIOLOGICAL EXAMINATION OF AN OBJECT 

DESCRIPTION 



This invention relates to a method for multiple-energy x-ray examination of an 

object. 

Radiological processes consist of having an x-ray pass through an object to be 
analyzed in order to deduce the distribution of various categories of materials in this 
object, absorbing the x-ray in different ways. A very common application is 
osteodensitometry, where the mass and density of bone tissue in a patient are analyzed, 
distinguishing these tissues from soft tissues. 

It is typical to use a broad spectrum of x-rays and to divide it into bands measured 
separately by respective measuring channels. Since the coefficients factors of absorption 
or attenuation of the x-ray by a given tissue category differ for each of the bands, the 
theoretical problem comes down to solving a system of equations, the number of which is 
the same as that of the measuring bands and each comprising two unknowns (the 
thicknesses and the masses of soft and bony tissue passed through). The problem 
becomes possible to solve after calibrating, achieved by having the x-ray pass through 
various standards provided with known thicknesses of materials with absorptive 
properties similar to those of the measuring material, in particular plexiglass and 
hydroxyapatite to simulate the soft tissue and the bony tissue. The mathematical 
parameters of a model 



linking the attenuation measurements to the thicknesses of the materials can than be 
calculated. j 

It should be noted that in actual fact, the human body contains three main 
categories of tissue: bony tissue, lean tissue and fatty tissue, but only two of them are 
generally taken into consideration due to the difficulty of distinguishing these three 
categories in the measurements, such that lean and fatty tissue are deliberately taken into 
consideration together. Other methods are then applied to distinguish their proportions in 
the soft tissue. 

The breadth of the spectrum makes it possible to have a must larger number of 
measuring bands than would be needed to obtain a solution, and to use all of them to 
obtain more precise results by making use of all of the absorption data obtained. In the 
article "Measurement of bone material using a multiple-energy x-ray absorptiometry," by 
J. Swanpalmer, R. Kullenberg, T. Hansson, Phys. Med. Biol, Vol. 43, 1997: pp. 379-387, 
where there are 23 measuring bands and where the three categories of tissue are taken 
into consideration, it is proposed to combine the measuring groups in threes in every way 
possible to obtain 1771 (23 x 22 x 21 / 2 x 3) systems of three equations with three 
unknowns that yield as many groups of results. A criterion of choice should then be 
applied. The authors advise choosing as the true result the one that is at the mean value, 
or the median value, for the most important parameter, which can be the bony mass 
passed through. 



Contrary to the above recommendation, the authors take into consideration the 
three categories of tissue. That does not call into question the validity of their method, 
provided the patient is exposed to a much stronger radiation intensity to reduce the 
uncertainties of the measurements to acceptable proportions. One could also simply 
modify their method and apply it to measurements on only two tissues, or more generally 
speaking, two categories of materials. However, this method is flawed in that it does not 
take into account the noise on the measurements, and it thus yields results with noise. 

The purpose of the invention is to perfect such methods by combining numerous 
results, and it comprises improving the result selection criteria. More precisely, in its 
most general form it relates to a method for radiological examination of an object in 
which at least two categories of materials are taken into consideration, comprising: the 
use of broad spectrum x-rays; measurements of the x-rays by bands of the spectrum; 
expressions (M?)of thicknesses or masses of the two categories of materials passed 
through by the x-rays, the expressions (M?) being functions of at least two of the 
measurements (mes k ) and coefficient (A); and applying a selection criterion from among 
the expressions (M?) to deduce from it an expression (final M?) considered true; 
characterized in that the selection criterion comprises a combination (f) of the expressions 
with weighting factors (a), and a 



calculation of weighting factors such that the combination has minimal noise (minimal 
variance in mathematical language) calculated according to the noise levels on the 
measurements (variance on the measurements). 

The invention will now be described referring to the figures: ■ 

- Fig. 1 is a spectrum view, 

- Fig. 2 is a distribution of results, 

- and Fig. 3 is a flowchart summarizing the invention, which can be referred to 
throughout the following description. 

The attenuation of the x-rays can be expressed as a function of the thicknesses 
passed through for each of the materials of indicators x and y, or of their masses M 
(density per surface unit) in the direction of the x-rays. The spectrum of measurements in 
Fig. 1 is divided into N bands generally marked by indicators i and j. The attenuations 
will vary in each of the bands according to variable absorption coefficients for the two 
materials. If we call mesj or mesj the measurements for a band of energy i or j, the masses 
passed through M x and M y can each also be expressed by the general formula 
(M?) = Ai + A 2 • mesi + A 3 • mesj + A4 • mes 2 i + A 5 • mes* + A 6 • mesj • mesj 
The measurements considered in this example being attenuation measurements, we will 
have for each measurement channel i (corresponding to a band of the spectrum ) the 
relation mesi = In (noi/ni) where noi is the number of photons arriving on the object and 
ni the number of photons having passed through the object. 



Since the non-linearity of the functions M x and M y as a function of the measurements is 
slight in practical reality, one can make use of this second degree polynomial function 
that comprises six coefficients Ai through A 6 . The degree of the polynomial can be 
adjusted depending on the problem. For example, to analyze objects made of materials 
having atomic numbers higher than that of biological tissue, as in non-destructive 
monitoring for examining a metal object, the measurements considered in this example 
being attenuation measurements, we will have for each measurement channel i 
(corresponding to a band of the spectrum ) the relation mesi = In (noi/ni) where noi is 
the number of photons arriving on the object and ni the number of photons having passed 
through the object 

These coefficients can be found in a phase of calibration through standards , 
sometimes called phantoms or wedges in the prior art, consisting of shaped parts with 
known, differing thicknesses made of materials simulating with their attenuation 
properties the materials of the object that will actually be measured. Each of these 
standards is thus radiated with the x-rays for a long period, making it possible to reduce 
the noise's influence on the measurements. The spectrum measured for each of the 
standards still yields N measurements resulting from the decomposition of the spectrum 
into as many bands. By then combining two series of measurements i and j taken for two 
bands N for each of the standards, the coefficients A are analyzed to adjust the functions 
M x and M y to the measurements. In this case, where there are six coefficients A for each 



of the two functions, and where the combinations of two bands of measurements are 
being considered, the measurements for each of the bands will have to tje for six 
standards in order to produce a unique solution. More numerous standards may also be 
used to improve precision on the solution. An error function minimizing will then be 
applied. 

This determining of coefficients A is repeated for various combinations of pairs of 
measurements. This was still the case in the above-mentioned previous article; however, 
it was noted that it was pointless to do all of the combinations, numbering Nx(N-l) 1 2, to 
fully utilize the measurements and that (N- 1 ) combinations were actually enough to 
gather all of the data. 

A preferred method consists of selling in the beginning the band of measurements 
with the least noise (for example the one that has the greatest signal n - number of 
photons at reception) and successively combining it with each of the other bands of 
measurements for the combinations. One finally obtains (N-l) estimates of the two 
functions M x and M y , which are noted M } ? , M 2 ? , M?N-\ for each of these two 
functions. 

At this stage of the process, the functions M x andM y representing the lengths 
passed through in two materials representing bony tissue and soft tissue can be converted 
to functions M u> M v and M w representing equivalent lengths passed through in the bony 
tissue, lean tissue and 



fatty tissue by combining linearly M x and M y in three different ways determined by 
experiment. This conversion method is independent of the invention anji already ready 
known in the prior art. 

Since there is no reason to prefer either of these M? estimates, a selection 
criterion must be applied to obtain the final M? estimate that will be considered true. In 
the previous article, one of the expressions obtained was directly selected according to a 
criterion of classification (the mean result) or median of the values considered by the 
expressions for one of the results. In the invention, the M? expressions will be combined; 
e.g. linearly according to the final M? formula = (al M?\) + (a2 M?2) +. . .+ (a N-l M7N 
-1) while minimizing the noise; the coefficients a is etc. have a sum equal to the unit 
(al +a2 + ... - + aN-l = 1). 

For each measurement channel, the noise on the number of photons follows a 
Poisson statistical law, the result of which is independent for each of the channels. The 
covariance matrix of the N-l results can be expressed according to the formula 

^N-l dM ? i DM ? j 1 

1\a = / x x ; 

^k=i 5mes k dmes k N K 

The variance on the linear combination yielding final M? is expressed by the 
formula f = (a i? a N -i) • T • l (ai, a N _i); this quantity f attains an optimal value when its 
derivative according to all of the variables 
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is null, that is, the noise's influence is minimized, i.e., 




af/aa N .! = o 



df/da 2 = 0 



The numerical resolution of this system provides the coefficients al, a2, etc. N 
final M? 9 that is, the masses passed through in the two categories of materials. 

These are the operations completed in the invention; it should be noted that the 
measurements (mes k ) used in the formula yielding r 5j and so forth are the measurements 
carried out through the actual object to be analyzed, but not those that were made by 
calibration to determine the coefficients A of the functions M?. 

The invention makes it possible to avoid the dilemma of prior methods in which 
part of the measurement energy was given up or, on the contrary, measurements that 
barely differed were accepted: it indeed uses the entire spectrum, but divides it into bands 
that are numerous enough to enable the measurement of each one to be usefully 
compared to other measurements made on distant bands. It thereby makes it possible, 
among other advantages, to analyze lean organisms as well as fatty ones. 

Lastly, the invention can be broadened to a greater number of materials than two, 
which can be advantageous in particular for imaging processes with contrast products, 
where three variables must be 



considered or for checking baggage (searching for explosives); and it can be applied 
when considering combinations of the M? functions other than linear orjes, 
Applications of the invention are: 

- osteodensitometry 

- obtaining the density of bone mass 

- obtaining body composition (fat mass, lean mass) 

- food-processing inspection, e.g. detecting bone fragments in meat or detecting 
pieces of glass in meals being served, 

- baggage inspection; searching for explosives, illegal products (weapons, food, 
drugs,...). 



